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SUMMARY 

The  definitive  orbit  for  Ariel  2  (l 964- 1 5A)  is  computed,  from  Minitraok 
observations,  for  a  period  of  twelve  months  from  the  launch  of  the  satellite. 
The  orbit  is  described  by  a  model  with  eight  orbital  parameters  and  these 
parameters  are  listed  at  every  twenty-fifth  nodal  passage.  The  angular 
observations  axe  accurate  to  about  1  *  and,  as  a  result,  the  average  computed 

standard  deviations  of  the  eight  fitted  orbital  parameters  are  as  follows: 
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1  m  in  semi-major  axis,  10  in  eccentricity,  2"  in  inclination,  4"  in  right 

ascension  of  the  node,  30"  in  argument  of  perigee,  0  ,03  in  time  at  the  node, 
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and  0,001  deg/d  and  0,001  deg/dJ  in  the  linear  and  quadratic  coefficients 
occurring  in  the  mean  motion  polynomial, 

Ephcmerides  computed  from  the  listed  orbital  parameters  will  bo  accurate 
to  about  -5-  km,  the  accuracy  required  by  the  Ariel  2  experimenters.  Limitations 
which  prevent  the  accuracy  from  being  better  than  this  arc  discussed. 
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1  INTRODUCTION 

Ariel  2  is  the  second  of  the  series  of  satellites  being  launched  in  the 
scientific  programme  based  on  Anglo-American  co-operation.  As  with  Ariel  1 , 
the  first  of  the  series,  the  satellite  was  constructed  and  launched  by  the 
U.S.  (NASA),  V/hereas  British  responsibility  for  Ariel  1  was  confined  to  the 
scientific  experiments  (including  telemetry  data  analysis),  for  Ariel  2  it  ha3 
been  extended  to  the  determination  of  definitive  orbital  parameters. 
Responsibility  of  the  U.K.  will  continue  to  grow  for  the  third  satellite  of  the 
series,  at  present  styled  U.K,  3  ani  due  to  be  launched  early  in  1967,  since 
this  will  be  the  first  spacecraft  actually  to  be  built  in  Britain, 

Three  experiments  constituted  the  scientific  payload  of  Ariel  2,  each 
sponsored  by  a  scientific  establishment  in  the  U.K. :  measurement  of  galactic 
noise  by  the  Mallard  Radio  Astronomy  Laboratory,  Cambridge;  measurement  of 
micrometeorite  flux  and  particle  sizes  by  the  Nuffield  Radio  Astronomy  laboratory, 
Jodrell  Bank;  and  measurement  of  atmospheric  ozone  by  the  Meteorological  Office, 
Bracknell.  The  satellite,  known  before  launch  as  S„52  or  U.K, 2  and  after  launch 
as  Ariel  2  or  1964-15A,  was  successfully  placed  in  orbit  at  17.25  UT  on 
27th  March  1964  from  tho  NASA  Wallops  Station,  Virginia,  by  a  four-stage  solid- 
propellant  rocket.  The  experiments  continued  to  work  until  the  end  of 
September  19^4,  by  which  time  the  spin  rate  of  the  satellite  wa3  too  low  for 
useful  data  to  be  obtainable. 

Tracking  data  have  been  provided  by  the  Mini  track  network  of  NASA  ( ST  AD  AN)  , 

The  first  observations  were  from  Lima,  Peru,  within  4  hours  of  the  time  of 

launch,  and  data  are  still  (November,  19^5)  being  obtained  from  tho  network, 

though  it  is  expected  that  NASA  will  officially  teiminate  the  project  in  the 

near  future.  The  data  have  been  analysod  in  Space  Department,  R.A.E. ,  and  the 

main  purpose  of  this  Report  is  to  tabulate  the  orbital  parameters  derived.  These 
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parameters  have  been  calculated  by  use  of  the  standard  programmes  ’  for  the 
Pegasus  computer  and  are  given  at  intervals  of  25  ascending  nodes  (about  I4  days). 
Computation  of  orbital  parameters  was  stopped  after  the  orbit  had  been  analysed 
for  twelve  months. 

2  MINI  TRACK  OBSERVATIONS 

The  STADAN  network  consists  of  a  dozen  Mini  track  stations  distributed  as 
shown  in  Pig.1.  The  inclination  of  the  Ariel  2  orbit  being  too  low  for  data  from 
Alaska,  observations  were  obtained  from  the  eleven  stations  listed  in  Table  1. 

The  Table  gives  latitude,  longitude  and  height  relative  to  the  Fischer 
ellipsoid**,  the  currently  best  world-wide  geodetic  datum.  Also  given  are  the 
nunfcer  of  observations  used  from  each  station. 
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A  total  of  about  3700  observations  has  been  used  over  the  twelve  month 
period,  corresponding  to  about  ten  per  day.  Of  these,  about  100  were  rejected 
during  analysis  due  to  the  size  of  residuals,  but  some  of  the  rejections  were 
due  to  incorrect  punching  of  capes  before  analysis  so  the  data  are  in  fact 
very  reliable. 

An  accuracy  of  the  order  of  a  millisecond  is  claimed  for  the  quoted 
times  of  observations,  and  time  errors  have  been  ignored  in  the  analysis. 

5 

Experience  with  Ariel  1  showed  that  the  angular  accuracy  is  effectively  about 
1  minute  of  arc  and  this  value  was  used  in  the  analysis.  It  covers  errors 
in  the  basic  interferometer  measurements,  ionospheric  refraction  correction, 
station  co-ordinates  and  orbital  model,  and  can  also  be  regarded  as 
incorporating  the  explicitly  ignored  time  errors. 

For  the  first  36  hours  after  launch,  observations  with  elevations  as  low 
as  20°  were  made  and  have  been  used  to  obtain  the  first  set  of  orbital  parameters. 
To  minimise  unknown  refraction  errors,  data  made  available  after  this  initial 
short  period  have  been  confined  to  observations  with  high  elevations;  more  than 
half  the  observations  have  had  elevations  greater  than  80°  and  there  have  been 
only  four  with  elevations  less  than  60°. 

The  coverage  of  each  interferometer  essentially  consists  of  a  pair  of 
narrow  fans  in  the  vertical  plane,  one  north-south  and  the  other  eaBt-west,  only 
one  being  used  on  a  given  pass.  For  th;Ls  reason  the  azimuths  of  most  observations 
have  been  close  to  0°,  90°,  180°  or  270°.  A  station  supplies,  for  a  given  pass, 
one  or  two  observations.  If  there  are  two,  since  they  are  both  in  the  same 
narrow  fan,  they  are  only  about  10  to  15  seconds  apart  in  time. 

3  ANALYSIS  OF  THE  OBSERVATIONS 

3.1  Dynamic  model  of  the  orbit 

The  analysis  has  been  carried  out  using  the  standard  R.A.E.  oibit 
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improvement  technique'  based  on  Merson's  smoothed  elements  .  The  orbital 

model  chosen  for  Ariel  2  has  twelve  parameters  associated  with  each  of  the 

epochs  defined  by  every  25th  ascending  node.  Eight  of  the  twelve  parameters 

are  determined  by  differential  correction  of  approximate  values;  these  eight  are: 

aQ  (semi-major  axis),  eQ  (ecdentricity),  iQ  (inclination),  0Q  (right  asoension 

of  the  defining  node),  u»o  (argument  of  perigee),  tQ  (time  at  the  defining  node) 

arel  the  parameters  n^  and  rig  such  that  the  mean  motion  at  time  t  is  given 

by:- 

n  =  (p/a^  +  nj  (t  -  tQ)  +  (t  -  to)2  , 
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where  p  =  398  602  km /see  .  The  remaining  four  parameters  -  denoted  by 
e,,  i. ,  Q.  and  w.  -  give  time-linear  contributions  to  e,  i.  fi  and  u;  unlike 
n^  and  n^,  they  are  not  fitted  to  the  observations  but  are  given  p  re -calculated 
values.  They  represent  drag,  luni-solar  and  earth-gravity  (Jy  and 
terms)  perturbations  of  the  orbit,  and  are  valid  for  only  two  or  three  days 
each  side  of  the  epoch  tQ.  It  must  be  anphasised  that  and  do  not  include 
the  large  secular  terms  which  arise  from  the  perturbation;  these  are  auto¬ 
matically  allowed  for  by  the  model.  Similarly,  e^  does  not  represej  t  the  main 
variation  of  e  due  to  drag,  this  being  allowed  for  through  the  parameter  n^. 
Some  printed  sets  of  parameters  have  been  circulated  in  which  coefficients 
02  and  (Og  were  included  in  addition  to  and  However  their  values  were 

so  small  as  to  be  quite  without  significance  (see  also  Section  5)» 

Ref. 3  gives  the  complete  set  of  formulae  by  means  of  which  geocentric 

position  is  computed  from  the  twelve  orbital  parameters.  They  differ  from  the 
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formulae  of  the  original  Pegasus  programme  in  two  respects.  First,  the 
position  of  the  satellite  at  time  t  is  computed  as  if  the  time  was  actually 
t  +  6t  where 


6t  =  0.088  sin  2  (0  -  0  -  18°)  seo  , 

0  being  the  sidereal  time.  This  formula  gives  a  first-order  approximation  to 
the  12-hour  along- track  oscillation  of  the  satellite  caused  by  the  tesseral 
harmonics  -  in  particular  by  J„  -  of  the  earth's  gravitational  field.  Further 
details  are  given  in  the  Appendix.  Naturally,  the  topocentric  position  of  the 
satellite  i3  found  by  taking  the  station  co-ordinates  (in  space  axes)  at  the  true 
time  t. 

The  other  divergence  from  the  formulae  of  the  original  Pegasus  programme 
is  that  a  term  has  been  introduced  into  the  expression  for  mean  anomaly,  M, 
associated  with  the  parameters  and  The  full  formula  for  M  is 

M  =  M0  +  (n0  -  U1  -  008  i0)(t  -  *0)  +  i  ^  (*  -  *o)2  +  i*>2  (t  -  to)3  , 

whore  U  is  determined  from  e  and  u  ;  t  is  taken,  of  course,  to  include  the 

O  O  07  '  ,  j_  ’ 

12-hour  6t  term  already  mentioned,  n  denotes  (p/sr)2  and  the  formula  for  n 
given  oarlier  is  not  affected. 

3.2  Computer  operation 

From  the  receipt  of  observations  from  NASA  to  the  generation  of  a  given 
set  of  definitive  orbital  parameters,  thore  are  seven  programmes  which  have  to 
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be  run  on  the  Pegasus  computer.  Allowing  for  the  fact  that  the  differential 
correction  programme  is  used  twice,  the  seven  programmes  correspond  to  eight 
operating  stages  as  listed  below,  The  first  two  stages  apply  to  the 
observations  in  general  and  the  remaining  six  to  the  analysis  at  a  particular 
node. 

(i)  The  direction  cosines,  of  which  the  data  consist,  are  punched  on 
paper  tape;  they  are  converted  to  azimuth  and  elevation  by  the  computer,  which 
provides  output  in  a  form  suitable  for  re-input  with  the  next  programme, 

(ii)  Observations,  sometimes  several  hundred  at  a  time,  are  stored  on 
the  Ariel  2  I.i, T. L.  (magnetic  tape  library). 

(iii)  The  initial  parameters  at  a  given  node  are  predicted  by  extrapolation 
from  the  definitive  parameters  25  nodes  before  (except  for  the  first  definitive 
node  when  NASA  estimates  are  used). 

(iv)  About  30  observations,  chosen  by  the  operator  for  the  analysis  at 
the  given  node,  are  transferred  from  the  M.T.L.  to  an  Ariel  2  S.D.T.  (selected 
data  tape);  some  preliminary  processing  is  incorporated  in  the  selection 
programme. 

(v)  A  single  iteration  of  the  main  differential  correction  programme  is 
run;  the  parameters  -  and  in  particular  n^  -  are  then  sufficiently  accurate  for 
the  next  stage, 

(vi)  Contributions  to  e^,  i^,  0^  and  due  to  gravity  ( J^,  and 
terms)  and  drag  perturbations  are  found;  these  are  functions  of  the  orbital 
parameters  and  the  drag  terms  depend  on  , 

(vii)  The  luni- solar  contributions  to  e^,  i_j,  0^  and  are  found  and 
added  in. 

(viii)  As  many  further  iterations  as  necessary  (normally  two)  of  the 

differential  correction  programme  are  run,  e^  etc  being  held  fixed;  after 

the  last  iteration  the  standard  deviations  of  a  ,  e  .  i  ,  H  ,  «  ,  t  ,  n.  and 

o'  o'  o'  o'  o'  o'  1 

ng  are  obtained. 

Stage  (ii)  requires  virtually  the  same  computer  time  -  about  10  min  - 
irrespective  of  the  number  of  observations  involved  (within  the  capacity  of 
the  M.T.L,).  For  the  other  stages  the  average  times  at  the  computer,  based 
on  a  set  of  30  observations,  are  as  follows:  for  stage  (i)  5  min;  for  (iii) 

2  min;  for  (iv)  8  min;  for  (v)  15  min;  for  (vi)  1  min;  for  (vii)  1  min;  for 
(viii)  30  min.  Thus  the  total  conputing  time  for  the  analysis  at  each  node 
worked  out  at  about  70  min  under  ideal  conditions.  An  actual  average  was 
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nearer  90  min  allowing  for  delays  due  to  the  rejection  of  observations  etc. 

The  time  away  from  the  machine,  including  a  little  more  than  1  hr  far  punching 
30  observations,  was  about  another  90  min.  Since  analysis  was  conducted  at 
210  different  nodes  the  routine  work  required  grand  totals  of  about  320  Pegasus 
hours  and  650  man  hours. 

3.3  Miscellaneous  points 

Times  of  observations  were  provided  relative  to  the  ’<7, TV  time  transmissions 
from  America.  Since  the  experimenters  were  using  telemetry  data  relative 
to  VAW  it  was  decided  not  to  correct  to  any  of  the  standard  time  systems  (universal, 
ephemeris  or  atomic  time).  If  it  is  so  desired,  the  orbital  parameters  in  this 
paper  may  be  corrected  to  datum  UT2  as  follows,  only  0q  and  tQ  being  affected: 
to  tQ  add  the  time  (in  seconds)  by  which  W..'V  is  slow  on  UT2;  to  Qq  add  0°.004* 
(UT2-WOT)  where  the  time  difference  is  in  seconds. 

The  analysis  at  each  defining  node  has  used  about  30  observations, 
spanning  a  period  of  3  to  4  days  centred  on  the  node.  Since  the  defining 
nodes  are  only  1^  days  apart  there  is  considerable  overlap,  about  half  the 
observations  being  used  twice.  This  was  done  deliberately  since  in  any 
application  of  the  'fitting  by  moving  arcs'  method  a  narrow  overlap  is  a  source 
of  inaccuracy.  Accuracies  actually  obtained  during  the  periods  of  overlap 
are  considered  in  Section  5. 

A  difficulty  arose  in  the  analysis  of  the  same  observations  at  two 
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different  nodes  since  WW  time  was  set  back  0  .1  at  1965  January  1.0.  Each 
observation  close  to  this  time  of  discontinuity  had  to  be  duplicated  in  the 
library  of  Ariel  2  observations,  once  with  the  original  observation  time 
and  again  with  an  adjusted  time,  the  observation  selected  depending  on  whether 
the  defining  node  used  was,  respectively,  the  same  side  of  the  discontinuity 
as  the  observation  or  the  opposite  side, 

4  RESULTS 

Orbital  parameters  for  Ariel  2  are  listed  in  Table  2.  Successive  columns 
of  the  Table  are  as  listed  below,  zero  suffixes  being  omitted  for  convenience. 

Node  number,  m 

Date 

Time 

Semi -major  axis,  a  (km) 

Eccentricity,  e 

10000  e^,  with  e^  in  units/100  days 

a  (1  -  e) 
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Inclination,  i  (degrees) 

100  i^,  with  i^  in  degrees/100  days 

Right  ascension  of  the  node,  fi  (degrees) 

100  0^,  with  01  in  degrees/lOO  days 

Argument  of  perigee,  u  (degrees) 

(degrees/lOO  days) 

Mean  anomaly,  M  (degrees) 

Mean  motion,  n  (degrees/lOO  days) 

n.  (degrees/(lOO  days)  ) 

1  —2  3 
10"  r^,  with  in  degrees/(lOO  days)"' 

Number  of  observations  used,  N 

Extent  of  the  observations,  D  (days) 

Standard  deviation  of  an  observation  of  unit  weight,  e 

Modified  Julian  Day  representation  of  date/time,  MJD* 

The  nodes  are  numbered  such  that  m  would  have  been  aero  at  a  fictitious 
node  half  an  hour  before  launch.  The  quantities  a(’i  -  e),  M  and  n  are  useful 
derived  parameters.  Although  a(l  -  e)  gives  an  indication  of  the  perigee 
distance,  r  ,  it  is  not  exactly  equal  to  r  •  To  0(J2)  the  perigee  distance  i3 
givjn  by:- 

2 

r  =  a  (l  -  e)  +  J  —  (sin2  i  cos  2  u  -  (l  -  e)(2  -  3  sin2  i)J  . 

The  values  of  the  eight  differentially-corrected  parameters  are  followed, 
in  Table  2,  by  their  computed  standard  deviations,  the  unit  in  each  case  being 
that  of  the  final  figure  quoted  for  the  parameter.  The  quantity  e  is  a  measure 
of  the  goodness  of  fit.  Its  expected  value  i3  1  (non-dimensional),  corresponding 
to  the  adopted  weighting  of  observations  with  standard  deviation  1 ' ;  if  e  is,  say, 
2.5  the  fit  is  such  that  the  observations  must,  a  posteriori,  be  regarded  as 
only  accurate  to  2^'.  Every  standard  deviation  includes  e  as  a  factor.  If  this 
factor  were  removed,  the  remaining  quantity  would  be  independent  of  fit  and 
would  depend  on,  and  therefore  provide  an  indication  of,  the  geometrical  coverage 
of  the  orbit  by  the  relevant  station  observations. 

It  is  seen  that  average  values  for  the  standard  deviations  of  the  eight 
basic  parameters  are  as  follows:- 
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a 

a 

1  metre 

°w 

o 

o 

• 

8 

~v) 

a 

e 

0.00001 

o 

o 

o 

V_M 

a. 

i 

0°.0005 

a 

n1 

12°/(100  days)2 

°h 

0°.001 

a 

"2 

Ip00°/(100  days) 

The  x’emarkable  accuracy  assessed  for  semi-major  axis  arises  from  the  relation 
n  a  =  |J  (Kepler's  third  lav/).  It  is  really  n  that  is  being  measured  so 
accurately  and  the  values  of  a  may  be  systematically  in  error  by  as  much  as 

x  2 

20  m  due  to  error  in  the  adopted  value  for  p,  viz,  398  602  km /sec  •  The 
above  c  corresponds,  in  fact,  to  (T  =  0°. l/lOO  days  or  a  maximum  contribution 
to  a„  (when  li|  days  av/ay  from  tQ)  of  about  0  .  002,  A  better  assessment  of 
orbital  distanoe  accuracy  is  given  by  cr  (a  1  -  e)  =  70  metres.  It  can 
readily  be  checked  that  the  values  given  for  and  cr^  also  both  correspond, 

for  t  -  t  =  I5  days,  to  contributions  of  about  0°.002  to  V  Thus  an 
approximate  estimate  of  the  maximum  along- track  error  in  position  computed  from 
the  definitive  orbital  parameters  may  be  obtained  from  S 3  x  0°.OO2  at  7000  km; 
this  estimate  is  les3  than  j  km.  further  discussion  of  errors  in  computed 
position  is  to  be  found  in  Section  5» 

To  indicate  the  behaviour  of  the  basic  parameters  relative  to  their 
tabulated  standard  deviations,  Pigs, 2-10  have  been  plotted.  The  idea  behind 
these  plots  has  been,  as  far  as  possible,  to  use  a  scale  large  enough  for 
the  representation  of  each  fitted  value  of  a  parameter  by  a  vertical  line  of 
length  2  a  centred  on  the  fitted  value.  To  this  end  secular  terns  have  been 
removed  from  some  of  the  parameters  by  preliminary  fitting  of  some  suitable 
polynomial,  and  for  two  of  the  parameters  -  e  and  w  -  the  long-periodic 
variation  lias  also  largely  been  removed. 

For  i,  n^  and  a  suitable  scale  is  available  without  modification  of 

the  data  from  Table  2;  Figs,5,  9  and  10  represent  these  three  parameters,  the 

plotted  lines  being  of  length  2  a. ,  2a  and  2a  respectively.  Eccentricity, 

1  n1  °2 

right  ascension  of  the  node  and  argument  of  perigee  have  been  represented  in 
Figs, 4,  6  and  7  by  graphs  of  Ae,  AQ  and  Aw  respectively,  where 
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Ae  =  e  +  1.08  x  10”^  m  -  7  x  10  ^  sin  u  , 

AO  =  n  +  k  x  360°  +  0°. 287515  m  +  2°.29  x  10“7  m2  +  1°.l8  x  10”11  m 

and 

Aw  =  (0  +  k  x  360°  -  0°,  21 385  m  -  2°,  3  x  10  7  m2  -  0°.4  cos  w  ; 

here  k  is  an  integer  which  has  been  introduced  to  avoid  the  discontinuities  of 
360°  in  the  tabulated  values  of  0  and  w. 

Figs. 2  and  3  both  relate  to  semi-major  axis.  The  values  of  <r  are  so 

cl 

small  that  the  scale  can  not  be  made  large  enough  over  the  whole  range  of  a, 
even  after  removal  of  a  very-high-order  polynomial.  Pig.  2  corresponds  to  the 
removal  of  a  quintic,  and  plots  A^a,  where 

A^  =  a  +  1. 61  x  ICf2  m  -  1.288  x  10"3  m2  +  6.915  x  10~9  m3  -  1.501  x  10"12  nA 

-16  5 

+  1.159  x  10  itT  . 

In  order  that  a  graph  showing  standard  deviations  might  be  drawn,  the  data  were 
divided  into  seven  sections  and  a  separate  polynomial  (quartic)  removed  from  each 
section,  the  constant  terms  being  adjusted  so  that  the  polynomials  joined  up. 

The  resulting  graph  i3  given  in  Pig. 3  by  the  plot  of  A^Q-,  where,  in  km, 

A^a  =  a  +  2.517  x  10  2  m  -  3.854  x  10  3  m2 

+  3.301  x  10"8  m3  -  1.082  x  lO"11  irA  for  0  <  m  <  750  , 

Aga  =  a  +  8.4985  -  9.22  x  10~3  m  +  1.539  x  10"3  m2 

-  6.976  x  10"9  m3  +  1.277  x  10~12  nA  for  750  <  m  <  1500  , 

A2a  =  a  -  180.9639  +  4.424  x  10"1  m  -  3.84  x  10”4  m2 

+  1.48277  x  10~7  m3  -  2.1107  x  10"11  nA  for  1500  <  m  <  2250  , 

A2a  =  a  -  1142.1030  +  1.83778  m  -  1.092788  x  10"3  m2 

+  2.88134  x  10’7  m3  -  2.82585  x  10”11  mA  for  2250  <  m  <  3000  , 
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A£a  =  a  +  43.1415  -  4.1146  x  ^0~2  m  +  2.0548  x  10-3  m2 

-  3.4637  x  10"9  m3  +  1.9394  x  10"13  m4  for  3000  <  m  <  3750  , 

=  a  +  375.2463  -  4.7361  x  10  ^  m  +  2.24885  x  10  4  ir2 

-  4.53198  x  10”8  ra3  +  3.34635  x  ICf12  m4  for  3750  <  m  <  4500  , 

and 

A2a  =  a  -  22497. 1421  +  18.545985  m  -  5.72321 6  X  10“3  m2 

+  7.845224  x  10"7  m3  -  4.027167  x  10”11  m4  for  4500  <  m  <  5250  . 

Pig, 8  relates  to  the  remaining  parameter,  tQ.  This  is  plotted  in  the 
Modified  Julian  Day  number  form,  a  quartic  polynomial  being  removed  from  the 
values  of  Table  2.  Again  standard  deviations  oannot  be  represented  and  this 
time  they  are  so  small  relative  to  the  scale  of  the  graph  -  the  average  cr 

b 

-7  0 

being  3.5  x  10  day  -  that  division  of  the  data  into  sections  would  not  help. 

The  quantity  plotted  in  fig. 8  is  At  ,  where,  in  days, 

AtQ  =  tQ  -  7.0318  x  10  ^  m  +  4*1  x108m2  +  4.06  x  10~12  m3  -  1.27  x  10  ^  m4  . 

Although  the  polynomial  (and  sinusoidal)  expressions  removed  from  the 
orbital  parameters  have  been  listed  above,  it  is  not  claimed  that  the  poly¬ 
nomials  represent  the  real  behaviour  of  the  parameters  or  that  the  shapes  of 
the  residual  plots  are  significant.  As  already  stated,  the  polynomials  have 
only  been  removed  in  order  that  the  magnitudes  of  the  standard  deviations  should 
be  visible,  on  the  graphs,  against  the  background  of  the  fluctuation  of  the 
parameters.  The  graphs  indicate  certain  facts  about  interpolation  that  are 
discussed  in  the  next  paragraph. 

Suppose  that  the  orbital  parameters  corresponding  to  a  certain  node  had 
not  been  available  and  that  they  we  re  estimated  by  a  suitable  form  of  interpolation 
in  the  adjacent  sets  of  parameters.  Then  Pigs, 4,  5,  6  and  7  show  that  the  values 
interpolated  for  e,  i,  Q  and  w,  respectively,  would  be  very  good;  it  would  be 
almost  certain,  in  fact,  that  any  interpolated  value  would  be  within  two  average 
standard  deviations  of  the  true  parameter.  For  a  and  n^,  however,  it  is  clear 
from  Pigs. 3  and  9  that  such  accurate  interpolation  would  be  impossible  and  for 
t  the  situation  would  be  worse  still.  (The  simplest  thing  to  do  with  the  final 
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parameter,  n^,  would  be  to  set  it  zero  -  Fig.  10  shows  that  the  error  would 
be  no  greater  than  in  interpolating  n^  •  ) 

Inadequacy  of  interpolation  for  the  parameters  a,  tQ  and  n«j,  due  to 
irregular  variation  in  air  drag,  explains  why  it  vra.s  necessary  to  obtain 
orbital  parameters  as  frequently  as  at  25  node  intervals  (see  also  Section  5) 
and  suggests  that  an  even  finer  analysis,  say  at  10  or  15  node  intervals, 
might  have  given  further  information  on  drag  fluctuations. 

5  ACCURACY  OF  POSITION  COMPUTATION 

Fran  a  set  of  orbital  parameters  derived  on  the  basis  of  a  given  orbital 
model,  the  position  of  a  satellite  at  any  time  may  be  computed  using  the  same 
orbital  model;  irrieed,  for  satellites  transmitting  scientific  measurements,  the 
provision  of  an  ephemeris  -  or  world  map  -  is  one  of  the  chief  motives  for 
determining  orbital  parameters.  In  the  case  of  .Ariel  2  the  experimenters 
wanted  positions  to  be  normally  accurate  to  j  km  or  better,  arid  almost  never 
worse  than  1  km. 

Successive  orbit  determinations  were  initially  carried  out  at  5G“no^® 
intervals,  starting  from  node  25.  Allowing  for  overlap  this  meant  that  a  given 
set  of  parameters  had  to  be  valid  for  a  period  stretching  from  about  26  nodes 
before  the  relevant  defining  node  to  about  26  nodes  after.  A  useful  measure 
of  the  accuracy  of  position  computation  was  available  by  comparing  the  positions 
at  various  times  -  in  particular  during  the  overlap  period  -  as  computed  from 
two  sets  of  parameters,  those  at  the  previous  defining  node  and  those  at  the 
following  one. 

Lot  d  be  the  distance  between  two  estimates  of  satellite  position  as 
above,  d  being  a  function  of  time.  The  behaviour  of  d  might  be  expected  to  be 
compounded  of  a  periodic  term  and  a  secular  term.  This  is  illustrated  in  Fig.  11, 
based  on  orbital  parameters  for  defining  nodes  75  and  125.  The  lower  plot  shows 
the  fine  behaviour  of  d  over  a  single  revolution  -  from  node  100  to  node  101. 

The  upper  plot  shows  the  coarse  behaviour  between  nodes  92  and  110;  a  curve 
has  not  been  drawn  because  of  the  oscillation  which  occurs  between  the  plotted 
points.  It  is  clear  at  once  that  the  desired  accuracy  i3  not  being  achieved. 

It  was  largely  as  a  result  of  contemplating  Pig,  11  that  the  decision  was 
made  to  produce  orbital  parameters  for  the  intermediate  nodes  50,  100  etc. 

Fig.  12  gives  the  coarse  and  fine  behaviour  of  d  for  compared  positions  based 
on  the  parameters  for  nodes  75  and  100.  Fig.  13  gives  a  similar  picture  based 
on  nodes  100  and  125.  Although  the  true  position  of  the  satellite  is,  of 
course,  never  k^own,  it  is  now  very  probable  that  the  desired  accuracy  is  being 
attained. 
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A  second  way  of  estimating  the  accuracy  of  position  computation  is  to 
employ  the  covariance  matrix  of  the  orbital  parameters  used  for  this  computation* 

A  Pegasus  programme  is  available  for  doing  this  and  has  been  described  elsewhere^. 
Let  S  be  the  root  mean  square  of  the  distance  between  the  true  position  of  a 
satellite  and  its  computed  position  from  the  parameters  at  a  defining  node  (only 
one  defining  node  is  involved  now).  Pig,  14  gives  the  coarse  behaviour  of  S 
between  nodes  1950  and  1963  based  on  the  covariance  matrix  of  the  orbital 
parameters  at  node  1950.  Pig.  15  gives  the  fine  behaviour  of  S  between  nodes 
1962  and  1963;  values  based  on  covariance  matrices  from  nodes  1950  and  1975  have 
been  averaged  here,  but  they  scarcely  differed.  On  the  same  graph  is  the 
corresponding  plot  of  d,  similar  to  that  plotted  in  earlier  figures. 

The  agreement  between  the  S  and  d  plots  of  Pig,  15  is  excellent,  bearing  in 
mind  that  they  arise  from  quite  different  methods  of  accuracy  assessment. 

Similar  plots  are  given  in  Pig. 16,  starting  from  a  different  pair  of  nodes,  2275 
and  2300.  This  time  the  agreement  is  much  worse,  an  average  value  for  d/S  being 
about  3|'«  If  the  errors  in  the  position  of  the  satellite,  computed  from  two 
sets  of  parameters,  were  uncorrelated,  the  expected  value  of  d/S  would  be  / 2, 
but  with  the  parameters  based  on  overlapping  sets  of  data  the  expectation  would 
be  less  than  this.  Thus  a  value  3'i  probably  represents  a  systematic  error  in  the 
orbited  model.  Hov/ever,  it  is  felt  tliat  even  in  this  case  the  accuracy  is  unlikely 
to  be  worse  than  1  1cm.  It  should  be  remarked  that  the  choice  of  nodes  upon  which 
Fig,  11  to  16  were  based  was  a  purely  random  one. 

The  effect  of  two  possible  sources  of  systematic  error  was  investigated  in 

a  similar  way  to  the  above,  using  the  quantity  d.  It  lias  been  remarked  in 

Section  3.1  tliat  parameters  0^  and  were  originally  included  in  the  model  but 

later  dropped.  The  effect  of  dropping  fl^  and  from  the  parameters  for  node  75* 

all  other  parameters  remaining  unchanged,  is  shown  in  Pig,  17.  This  effect  is, 

of  course,  proportional  to  (t  -  t  )  and  its  maximum  value  may  be  seen  to  be, 

2  ® 

in  metres,  about  33  (t  -  t  )  ,  measuring  time  in  days.  This  is  completely 
negligible. 

The  Fischer  ellipsoid  was  adopted  early  in  the  analysis,  but  without 
recomputing  all  previous  sets  of  parameters  obtained  when  the  stations  were  not 
all  referred  to  the  same  datum.  One  set  of  parameters  was  recomputed,  the  set 
for  node  525.  Pig.  18  indicates  the  negligible  change  in  position  as  computed 
from  the  parameters  before  and  after  the  change. 
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6  COMPARISON  iTITH  OTHER  ORBITS  OF  ARIEL  2 

Elements  of  Ariel  2  have  been  issued  by  two  American  computing  centres, 
Spacetrack  and  NASA*  In  each  case  new  sets  of  elements  are  produced,  on  average, 
about  every  nine  days,  some  five  times  less  frequently  than  the  elements 
tabulated  in  this  paper.  For  this  reason,  among  others,  the  American  orbits 
must  be  regarded  as  significantly  less  accurate  than  the  definitive  R.A.3. 
orbit.  The  Spacetrack  bulletins,  in  any  case,  describe  their  data  as  for 
prediction  purposes  only  and  unsuitable  for  scientific  analysis. 

The  NASA  elements,  like  the  R.A.S.  elements,  have  been  derived  from  Mini¬ 
track  observations.  However,  in  addition  to  being  issued  less  frequently  they 
differ  from  R.A.E.  elements  in  that  they  are  related  to  a  different  dynamic 
model  of  the  orbit;  also,  each  set  corresponds  to  some  arbitrary  epoch  and  no+' 
to  the  time  of  a  particular  nodal  passage.  The  precise  definitions  of  NASA 
elements  are  rather  difficult  to  discover  and  it  is  not  certain  that  the  same 
orbital  model  is  used  for  the  analysis  of  all  satellites.  It  was  therefore 
of  interest  to  compare  NASA  elements  -  interpolated  to  the  time  of  the  nearest 
set  of  R.A.E.  elements  -  with  the  R.A.E.  elements,  using  the  differences  to 
establish  the  most  likely  interpretation  of  the  NASA  elements. 

The  comparison  has  been  made  and  it  seems  that  NASA  must  have  used  a 
set  of  elements  which  will  be  defined  in  the  next  paragraph.  They  are  denoted 
by  a,  e  etc.  to  distinguish  them  from  R.A.E,  elements.  Formulae  are  given  in 
Sections  6.1  and  6.2  which  relate  the  two  sets  of  elements.  On  the  assumption 
that  the  NASA  elements  are  the  barred  elements,  a  number  of  sots  have  been 
transformed  into  Merson*s  elements  for  comparison  with  R.A.E.  elements.  The 
averago  difference  is  then  about  three  or  four  times  the  R.A.E,  standard 
deviations  and  -  since  the  R.A.E.  elements  are  the  more  accurate  -  the  agreement 
is  very  satisfactory.  This  is  illustrated  by  Fig.  19  which  gives  the  R.A.E. 
eccentricity  graph;  values  of  e  from  NASA  are  plotted  before  and  after  trans¬ 
formation  and  it  is  seen  that  after  transformation  they  lie  on  or  near  the 
R.A.E,  curve.  Comparison  of  R.A.E,  and  NASA  epheraorides  is  considered  in 
Section  6,3  and  illustrated  by  Fig,  20, 

The  elements  a,  e  etc.  are  defined,  to  0(e),  by  a  double  averaging  of 
osculating  elements.  The  first  averaging  is  with  respect  to  mean  anomaly  and 
the  second  with  respect  to  argument  of  perigee;  these  averagings  remove, 
respectively,  short-periodic  perturbations  associated  with  the  J2  term  of  the 
geopotential  function  and  long-periodic  perturbations  associated  with  the 
term.  Despite  the  general  uncertainty  over  their  orbital  methods,  there 
la  reason  to  believe  that  NASA  remove  perturbations  also,  but  these  are 
much  smaller  than  the  perturbations  and  are  not  considered  hero. 
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The  formulae  given  below  have  been  derived  by  using  the  results  of  Merson  , 

6  9 

who  compared  his  smoothed  elements  with  the  mean  elements  of  Kozai  *  The 
elements  e,  i  and  fi  are  very  close  to  those  which  Kozai  reached  after  removing 
long-periodic  terms  but  the  a  and  u  defined  here  differ  from  those  of  Kozai 
since  each  of  the  latter  was  given  a  bias  relative  to  the  true  average  short- 
periodic  perturbation. 

The  element  a  is  dealt  with  after  e,  i,  Q  and  u  in  the  comparisons  below 
and  is  then  considered  in  conjuration  with  the  mean  motion  n  and  the  anomalistic 
period.  Kean  anomaly,  for  which  comparison  would  not  be  so  straightforward,  is 
omitted  from  consideration, 

6. 1  Eccentricity,  inclination,  right  ascension  of  the  node  and  argument  of 
perigee 

To  0(e)  the  relations  between  e,  i,  Q,  u>  and  e,  i,  Q,  u  respectively  are 
given  by:- 

e  -  e  = 

i  -  I  = 

n  -  Q  = 

and 

0)  -  u  = 


The  first  term  in  the  expression  for  c  -  e  is  negligible  here  since 
2 

2-3  sin  i  =  0(e);  the  second  term  represents  the  transformation  made  to  the 
plot  of  NASA  eccentricity  in  Pig,  19. 

6.2  Somi-ma.ior  axis,  mean  motion  and  anomalistic  period 

The  semi -major  axis  and  the  mean  motion  arc  directly  related  by  Kepler's 
third  law.  Since  -  with  observations  of  direction  and  not  range  -  it  is  n  rather 
than  a  which  is  measured,  it  is  preferable  to  give  the  formula  relating  n  and  n. 
There  being  no  long-periodic  term,  the  relation  -  to  0(e)  -  is:- 


i 


/V\2  o  a  r 

J2  (r)  e  (2  "  3  sin"  *•)  -  2  jJT  sin  1  sin  “  * 


1  /p\2  1  *^5  ^ 

-  J2  ( — )  sin  2  i  (3  +  A  e  cos  to)  +  —  e  cos  i  sin  u  , 


2  J2  a 


1  AA2  1  ® 

—  <4,  f  rl  cos  i  (3  u  +  3  M  +  4  e  sin  u)  -  —  ■“ —  e  cot  i  cos  w 

O  ^ 


2  °2  a, 


2 

”  4  ^2  ^(4  -  5  sin2  i)(3  u  +  3  Mq  V  2  e  sin  w)  -  ^  sin2  i  sin  2  wj 


1  il  -i  .  . 

rr  t  —  c  sin  i  cos  u 

2  J2  a 
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n  -  n 


n  (2  -  3 


•  2  .s 

sin  i) 


Although  the  true  anomalistic  period  -  the  time  from  one  passage  through 
perigee  to  the  next  -  is  2fl/n,  it  appears  to  he  2ft/n  which  is  listed  as 
•anomalistic  period’  hy  NASA.  Thus  the  latter  is  consistently  half  a  second 
greater  than  the  true  anomalistic  period. 

Values  listed  in  the  Spacetrack  bulletins  do  appear  to  be  of  the  true 
anomalistic  period  and  agree  satisfactorily  with  values  of  2%/n  using  Table  2. 

6.  3  Comparison  of  ephemerides 

It  does  not  matter  that  the  definitions  of  R.A.E.  and  NASA  elements  are 
different,  so  long  as  the  appropriate  formulae  are  used  for  generating  satellite 
positions  in  each  case.  A  comparison  has  been  made  of  satellite  position 
computed  from  R.A.E.  orbital  elements  with  the  ephemoris  provided  by  the  NASA 
•Refined  World  Maps'.  The  period  from  node  600  to  node  601  was  chosen  for  the 
comparison  since  during  this  period  the  difference  between  e  (NASA)  and  e  (R.  A.E. ) 
was  near  its  maximum  (0.0000),  as  may  be  seen  from  Pig.  19. 

Fig.20  gives  (lower  curve)  the  difference  between  the  R.A.E.  and  NASA 
computed  satellite  positions  ard  (upper  curve)  the  height  component  of  this 
difference.  It  is  seen  that  the  height  difference  does  not  exceed  310  metres 
nor  the  total  difference  800  metres.  Since  the  R.A.E.  elements  are  the  more 
accurate,  position  agreement  to  this  order  is  very  satisfactory  and  confirms  that 
both  sets  of  orbital  elements  are  being  used  correctly. 

7  DISCUSSION 

The  accuracy  required  by  the  Ariel  2  experimenters  has  been  achieved  in 
the  orbit  provided.  Had  an  order  better  accuracy  been  required,  a  careful  analysis 
of  the  possible  sources  of  error  would  have  been  necessary.  At  the  expense  of 
some  complication,  further  improvements  in  the  dynamic  model  of  the  orbit  might 
have  been  mado  and  the  Pegasus  computer  programme  improved.  However,  it  is  doubt¬ 
ful  whether  the  Minitrack  data  themselves  would  have  been  accurate  enough  to  warrant 
much  refinement  in  the  computing. 

The  average  size  of  the  residuals  in  angle  is  the  1  minute  of  arc  that  had 
boon  expected  when  the  analysis  was  started.  For  a  satellite  at  an  average 
distance  of  about  1000  km  this  corresponds  to  an  error  in  distance  of  about  300  m. 

For  observations  with  elovations  as  low  as  60°  the  greatest  source  of  error 
is  probably  the  inadequate  correction  of  ionospheric  refraction.  It  has  been 
estimated*1^  that  for  an  elevation  angle  of  50°  the  residual  error  after  refraction 
correction  is  2  minutes  of  arc  so  that  for  60°  it  is  still  important. 
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Other  errors  in  the  data  are  not  likely  to  be  significant.  Systematic 
errors  which  arise  during  the  periods  between  successive  calibrations  of  a 
llini track  station  -  periods  of  some  months  in  general  -  should  not  average  more 
than  4  minute  of  arc  ,  Errors  in  station  survey  should  not  exceed  100  m, 
leading  to  errors  of  the  same  order  as  calibration  errors. 

Turning  to  the  orbital  model,  appreciable  errors  -  say  about  200  m  -  still 
arise  from  inadequate  allowance  for  the  earth's  tesseral  harmonics  (see  the 
Appendix),  despite  the  empirical  modification  made.  Errors  from  neglect  of 
the  zonal  harmonics  beyond  are  thought  to  be  negligible.  The  effect  of  Jy, 

for  example,  is  not  important  for  a  satellite  of  inclination  51°*6.  The  values 
of  e^  in  Table  2,  which  are  from  theory,  arc  in  general  agreement  with  values 
from  differencing  e,  bearing  in  mind  that  the  orbital  model  automatically  allows 
for  the  main  drag  term  in  e  -  i.e#>  if  is  zero  it  is  a(l  -  e)  and  not  e 

which  is  held  constant. 

Errors  from  inadequate  representation  of  the  effects  of  atmospheric  drag 
are  less  easy  to  estimate.  They  may  be  divided  into  short-periodic  and  long- 
periodic  errors.  The  short-periodic  errors  arise  because  drag  is  concentrated 
at  perigee  and  not,  as  the  mean  anomaly  polynomial  would  suggest,  uniformly 
spread  over  the  orbit.  But  the  errors  are  negligible,  being  less  than  a  second 
of  arc  -  in  the  case  of  Ariel  2  -  for  mean  anomaly  and  so  equivalent  to  at  most 
a  few  seconds  of  arc  in  an  observation. 

The  long-periodic  errors  in  the  representation  of  drag  effects  are  more 

serious  and  it  is  these  which  are  difficult  to  estimate.  They  arise  because  of 

fluctuations  in  the  density  of  the  atmosphere  which  are,  on  the  different  time 

scale  involved,  of  short  period.  These  fluctuations  have  been  investigated  by 

11 

King-Hele  and  Quinn  who  have  plotted  values  of  n^,  as  in  Pig. 9,  as  far  as 
10th  September,  1964*  Geomagnetic  index,  when  plotted  for  comparison  with  the 
parameter  n^ ,  shows  the  same  general  behaviour  but  has  very  pronounced  'spikes'. 
It  is  clear,  in  fact,  that  the  calculation  of  n^,  even  at  intervals  as  fine  as 
1a  days,  has  smoothed  the  true  orbital  acceleration  and  thereby  introduced  error. 
If  the  required  corrections  to  n^  could  be  estimated  and  doubly  integrated, 
then  the  correction  to  mean  anomaly  could  be  found  as  a  function  of  time.  Prom 
a  rough  inspection  of  the  n^  and  geomagnetic  index  curves  it  seems  likely  that 
the  error  could  reach,  or  exceed,  a  minute  of  arc  in  equivalent  topocentric 
observation. 

It  is  interesting,  in  the  light  of  the  last  paragraph,  to  speculate  that 
to  obtain  more  accurate  orbits  for  satellites  like  Ariel  2  it  may  be  necessary 
to  feed  in  values  of  geomagnetic  index  over  the  t hree-or-f our-day  period  of  an 
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orbit  determination  and  to  form  terms  contributing  to  the  mean  anomaly  by  a 
double  numerical  integration  of  these  values  of  the  index. 

In  addition  to  its  application  to  the  study  of  short-periodic  fluctua¬ 
tions  in  atmospheric  density  the  orbit  of  Ariel  2  lias  been  used  in  the 

evaluation  of  the  odd  zonal  harmonics  of  the  earth's  gravitational  potential. 

12 

King-Hele,  Cook  and  Scott  have  taken  Ariel  2  as  one  of  six  satellites  for 
this  evaluation.  They  used  the  data  for  perigee  radius,  a  (l  -  e),  rather 
than  eccentricity,  in  order  to  remove  drag  effects  as  accurately  as  possible, 
and  found  that  the  corrected  data  we  re  then  as  good  as  data  for  tho  other 
satellites  in  virtually  drag-free  orbits. 

8  CONCLUSIONS 

By  use  of  an  orbital  model  with  eight  parameters  fitted  'ran  Minitraok 
observations,  the  orbit  of  Ariel  2  (l  9&4-15A)  has  been  successfully  computed 
for  a  period  of  a  year  from  its  launoh.  The  orbital  parameters  have  been 
listed  (in  Table  2)  for  the  time  of  every  twenty-fifth  passage  through  the 
ascending  node. 

Ephemerides  calculated  from  the  listed  parameters  should,  in  general, 
bo  accurate  to  about  j  km,  i.e.  to  the  accuracy  required  by  the  experimenters. 
This  estimate  of  accuracy  is  based  on  the  view  that  the  effective  average 
error  in  a  Minitrack  observation  is  1  minute  of  arc,  a  figure  supported  by  the 
residuals  in  the  observations  after  orbits  have  been  fitted.  The  estimate  has 
been  confirmed  by  the  comparison,  over  a  short  period,  of  ephemerides  calculated 
from  two  sets  of  orbital  parameters,  corresponding  to  the  epochs  before  and 
after  tho  period  in  question. 

Consideration  has  been  given  to  the  factors  limiting  the  accuracy  of 
the  orbital  determinations.  It  is  concluded  that  the  throe  main  factors  are: 
inadequacies  in  the  correction  of  observations,  other  than  at  very  high 
elevation,  for  ionospheric  refraction;  inadequacies  in  the  representation  of 
tesseral  terms  in  the  earth's  gravitational  field;  and  inadequacies  in  the 
representation  of  drag  effects  over  periods  of  several  days, 
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Appendix 

PERTURBATIONS  OP  AREL  2  DUE  TO  THE  TSSSERAL 
HAR^OUICS  OP  TUB  BARTHS  GRAVITATIONAL  FIELD 

During  the  first  few  orbital  determinations,  based  on  the  standard  dynamio 

model  that  was  being  used  at  the  beginning  of  1964,  it  was  noticed  that  residuals 

from  the  V/inkf  ield  Minitrack  station  were  consistently  larger  than  residuals  from 

the  other  stations.  It  was  discovered  moreover  -  using  the  original  form  of  the 

orbit  determination  programme  in  which  residuals  in  observation  times  as  well 

2 

as  angular  data  were  computed  -  that  the  residuals  were  almost  entirely  along- 
track,  equivalent  to  an  error  in  the  V/inkf  ield  clock.  This  error  was  oscillatory, 
with  amplitude  about  0S.  1  and  period  just  under  (by  about  10  minutes)  half  a 
day.  It  was  realised  that  an  error  of  this  type  could  be  caused  by  the  neglect 
of  the  coefficient  (J  _)  of  one  of  the  tesseral  harmonics  of  the  earth* s 
gravitational  field,  though  it  had  earlier  been  thought  that  ^  Perturbations 
would  not  be  detectable.  The  proof  that  the  effect  was  in  fact  caused  by 
gravitational  perturbations  -  and  that  the  V/inkf  ield  clock  was  blameless  -  was 
obtained  by  plotting  time  residuals  for  all  the  185  observations  which  were 
obtained  from  the  Minitrack  network  during  the  first  90  hours  from  launch.  The 
residuals  were  plotted  against  the  argument  0-0,  where  0  is  sidereal  time  and 
0  the  right  ascension  of  the  satellite  node. 

Let  the  negative  of  the  residual  in  an  observation  time  (assuming  the 
original  form  of  the  orbit  determination  programme)  be  denoted  by  6t.  Then 
6t  is  the  amount  by  which  an  observation  time  has  to  be  increased  before  computing 
satellite  position  from  the  standard  dynamic  model.  It  was  found  that  a  good 
fit  to  the  plot  of  185  St's  was  given  by:- 


6t  =  0,088  3in  2(0  —  Q  -  18°)  sec  . 


Now  the  formula  for  the  along- track  effect  of  ^  given^  ^ >  for  a  circular 
orbit,  by:- 


6L  *  Had 


2  .  2 
n  sin 


sin  2(6  -  n  +  J 


2,2 


where  p,  n  and  i  are  semi-latus  rectum,  mean  motion  and  inclination  respectively, 
R  is  the  (mean)  equatorial  radius  of  the  earth,  \  _  is  the  (east)  longitude 
of  one  extremity  of  the  major  axis  of  the  earth's  equator,  assumed  elliptical, 
arxl  l  is  8  -  ft,  Taking  6L  =  n  6t  and  substituting  the  values  of  p,  n,  i  and  ft 
for  the  Ariel  2  orbit,  the  observed  6t  fit  may  be  accounted  for  by  a  J.  . 
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Appendix 


perturbation  for  which  J,  .  =  3.0  x  10  and  _  =  -18  .  The  orbit  deter- 
mination  programme  has  consequently  been  modified  by  arranging  that  the  time 
of  an  observation  is  increased  by  fit  where 


5t 


2 

=  2  J2,2 


.  2  . 

— “  sin  2(0-O  +  X22) 


t 


J_  „  and  \0  0  having  the  values  above. 

14- 

Now  the  values  of  J0  0  and  X0  Q  obtained  by  Guier  and  Newton  ,  with 

-6  o 

which  other  values  are  in  fair  agreement,  are  1.72  x  10  and  -13  *4  respectively, 
so  that  at  first  sight  the  correction  is  almost  double  what  can  be  justified. 

The  natural  explanation  is  that  the  fitted  J2  2  accounts  not  only  for  the  pure 
J2  2  but  also  for  (the  along-track  effects  of}  2,  Jg  2,  Jg  2  etc.  However, 
this  explanation  has  not  been  satisfactorily  borne  out  on  evaluating  these 
effects  a3  far  as  Jg  2,  using  values  from  Ref.  14.  ftie  pure  2  effect  and 
the  J^  2,  Jg  2  and  Jg  2  effects  are  then  given,  respectively,  byr¬ 


and 


fit  =  0.051  sin  2  (8  -  0  -  13°)  , 

6t  =  0.017  sin  2  (0  -  0  -  67°)  , 

6t  =  0.004  sin  2  (0  -  0  +  22°) 

fit  =  0.002  ain  2.(0  -  n  -  11°)  ; 


the  formulae  for  the  quantities  L_  and  L 

^,2  4|2 

similar  formulae  developed  for  Lg  2  and  Lg  2» 
given  by 


of  Ref,  13  have  been  used  here,  and 
On  combining,  the  total  effect  is 


fit  =  0.051  sin  2  (  6  -  0  -  21°) 


so  that  the  amplitude  is  the  same  as  for  the  pure  J2  2  effect  though  the  phase 
has  ohangod.  Values  for  J  _  „  etc.  are  not  available  but  there  seems  to  be 

I  £ 

little  hope  of  increasing  the  amplitude  to  the  observed  0.088  sec  through  their 
91 

agency. 


The  effioaoy  of  the  empirical  programme  modification  may  be  appreciated 
when  it  is  stated  that  the  mean  value  of  the  observation  residuals  is  halved. 
Although  the  modification  was  based  on  the  analysis  at  the  early  nodes  it  was 
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tested  also  at  later  nodes;  at  node  3250,  for  example,  the  mean  (rms)  residual 
was  2**80  without  the  modification,  but  1**46  with  it* 

It  is  recalled  that  the  programme  modification  accounts  for  only  that 
component  of  the  pertubation  due  to  tesseral  harmonics  which  manifests  itself 
as  an  apparent  along-track  error  in  satellite  position  with  period  about  half 
a  day.  There  are  also  components  in  the  two  directions  perpendicular  to  the 
along-track  direction.  The  perturbation  which  is  cross-track  but  within  the 
orbital  plane  contains  eccentricity  as  a  factor  and  may  be  neglected,  but  the 
perturbation  perpendicular  to  the  orbital  plane  is  significant.  Its  effect  is 
less  noticeable  than  that  of  the  along-track  component,  partly  because  it  has  a 
short-periodic  factor  superimposed  on  the  half-day  period.  Its  removal,  however, 
might  be  expected  to  contribute  to  a  further  improvement  in  orbital  fits. 

For  a  complete  representation  of  the  effect  of  any  particular  tesseral 
harmonic  the  natural  procedure  would  be  to  introduce  the  perturbations  of  the 
six  orbital  elements.  With  the  normal  notation  J  (but  with  <T  for  the  element 
sometimes  denoted  by  \*)  there  would  be  no  perturbation  in  a,  while  perturbations 
in  w  and  could,  for  a  near-circular  orbit,  be  combined;  thus  expressions 
for  6e,  5 i,  6Q  ard  6u  +  6c  would  be  required  instead  of  merely  the  expression 
for  5L  (=  6u  +  6c  +  60  cos  i). 


Attempts  have  been  made  to  represent  the  complete  perturbations  of  Ariel  2 
due  to  J„  »  and  J  „  by  taking  either  '"he  values  of  Ref,  14  or  else  the  empirical 

values  for  the  along-track  perturbation  only,  viz,  J2  2  =  x  10  * 

o  ^ 

„  =  —18  and  J.  „  =  0,  Both  attempts  were  unsuccessful  -  residuals  became 
4,4 

larger  not  smaller  -  ard  perhaps  this  is  not  surprising,  Since  it  was  necessary 
to  take  an  empirical  J ^  based  on  residual  fitting,  to  represent  the  along- 
track  perturbation  alone  successfully,  it  would  presumably  be  necessary  tc  take 
a  new  empirical  value  to  represent  the  complete  perturbation.  But  die  method 
of  residual  fitting  would  not  then  be  so  obvious;  since  the  revised  orbital 
model  -  with  the  empirical  along-track  correction  -  gave  the  Ariel  2  orbit  to 
the  required  accuracy,  the  question  of  further  improvement  has  not  been  closely 
studied* 


The  perturbations  due  to  J  with  s  *  2  would,  of  course,  have  to  be 

n,s 

considered  in  a  full  analysis.  The  tables  below  give  order-of -magnitude  estimates 
of  'fie  amplitudes  of  the  perturbations  from  ^  as  far  as  ^  using  values 
from  Ref.  14,  the  effects  on  the  elements  being  represented  in  metres.  Lines  of 
the  tables  for  which  the  perturbations  include  eccentricity  (e  =  0,07)  as  a 
factor  are  indicated  by  an  asterisk. 
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!  H  J  *  -UL'M  3a  Wsl  I  ilWfl  ■*? 


Station  location 

Blossom  Point,  Maryland,  USA 
Cast  Grand  Forks,  Minnesota,  USA 
Fort  Myers,  Florida,  USA 
Goldstone  Lake,  California,  USA 
Hartebeeshcek,  Johannesburg,  SA 
Lima,  Peru 
Quito,  Ecuador 
Santiago,  Chile 

i 

'  St.  Johns,  Newfoundland 
1  Ulnkfleld,  England 
Woomera,  Australia 


Latitude 

Longitude 

Height  (metres) 

No.  of  observations 

38,63053  N 

77.08629  W 

-1 

252 

68,02256  N 

97.01082  W 

256 

760 

26.5/1827  N 

81 .86539  w 

7 

225 

35.33016  N 

116.89977  W 

937 

202 

25.88361  S 

27.70791  E 

1571 

216 

11.77535  s 

77.15026  w 

2 

210 

oo.6>^07  s 

78.57900  w 

3568 

160 

33.16896  s 

70.66865  w 

636 

176 

67.76137  N 

52.72036  w 

106 

698 

51.66595  N 

0.69623  w 

91 

669 

31.39167  S 

136.86972  E 

118 

189 

Latitude  and  longitude  are  referenced  to  the 
Fischer  ellipsoid!  a  spheroid  of  semi-major 
ails  6378.166  km  and  flattening  1/298.3, 
Heights  are  measured  upwards  from  this 
spheroid. 
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ORBITAL  PARAMETERS  OP  ARIEL  2 
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5175  IUH  25  0  41  52.80  3  7150.8553  e  0.068095  a  -42  6663.92  51.6426  5  1  174.2194  a  26  170.551  7  -1.54  i-169.203  516855.48  2781  8  -83  m  56  3.3  1.4  38844.0290833 

5200  IUH  26  18  27  30.66  3  7150.4230  9  0.067954  6  -42  6664.52  51.6424  5  3  166.9447  9  25  175.947  5  -1.13  i -175. 367  516902.35  2587  9  -103  is  41  3.4  1.5  38845.7691048 
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5250  IUH  30  5  58  8.69  3  7149.695012  0.067697  a  -43  6665.68  51.6433  6  -2  1  52.3919  12  22  1  86.780  i  -0.45  1-187.744  516981.30  1  968  26  -91  39  31  2.8  1,3  38849.2487117 
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FIG.  8  EPOCH,  t0,  WITH  QUARTIC  POLYNOMIAL  REMOVED 
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FIG.  II  DIFFERENCE  BETWEEN  COMPUTED  SATELLITE  POSITIONS 
BASED  ON  ORBITAL  PARAMETERS  FOR  NODES  75  AND  125 
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FIG.  12  DIFFERENCE  BETWEEN  COMPUTED  SATELLITE  POSITIONS 
BASED  ON  ORBITAL  PARAMETERS  FOR  NODES  75  AND  IOO 
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FIG. 1 3  DIFFERENCE  BETWEEN  COMPUTED  SATELLITE  POSITIONS 
BASED  ON  ORBITAL  PARAMETERS  FOR  NODES  IOO  AND  125 
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FIG.  14  PLOT  OFS  BASED  ON  COVARIANCE  MATRIX  FOR  NODE  1950 
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FIG.  1 5  ACCURACY  OF  COMPUTED  SATELLITE  POSITION 


ASSOCIATED  WITH  PARAMETERS  FOR  NODES  1950  &  1975 
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FIG  16  ACCURACY  OF  COMPUTED  SATELLITE 
POSITION  ASSOCIATED  WITH  PARAMETERS 
FOR  NODES  2275  AND  2300 
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BOTH  PLOTS  BASED  ON  ORBITAL  PARAMETERS  FOR  NODE  52.5 


FIG.  1 8  VARIATION  OF  COMPUTED  SATELLITE  POSITION 
DUE  TO  ADOPTION  OF  FISCHER*  ELLIPSOID 
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FIG.  20  COMPARISON  OF  NASA  AND  RA£.  EPHEMERIDES 
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